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We derive the efTective Hamiltonian for a quantum system constrained to a submanifold (the 
constraint manifold) of configuration space (the ambient space) in the asymptotic limit where the 
restoring forces tend to infinity. In contrast to earlier works we consider at the same time the effects 
of variations in the constraining potential and the effects of interior and exterior geometry which 
appear at different energy scales and thus provide, for the first time, a complete picture ranging 
over all interesting energy scales. We show that the leading order contribution to the effective 
Hamiltonian is the adiabatic potential given by an eigenvalue of the confining potential well-known 
in the context of adiabatic quantum wave guides. At next to leading order we see effects from 
the variation of the normal eigenfunctions in form of a Berry connection. We apply our results to 
quantum wave guides and provide an example for the occurrence of a topological phase due to the 
geometry of a quantum wave circuit, i.e. a closed quantum wave guide. 
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I. INTRODUCTION 

The derivation of effective Hamiltonians for constrained 
quantum systems has been considered many times in the 
literature with different motivations and applications in 
mind. Roughly speaking, the available results split into 
two different categories which are related to two different 
energy scales. In the context of adiabatic quantum wave 
guides one considers the situation where the strong forces 
restricting the particle to the wave guide change their 
form along the direction of propagation. The eigenvalues 
of the transverse Hamiltonian thus also vary along this 
direction and produce an effective adiabatic potential for 
the tangential dynamics, i.e. for the propagation. In this 
case the tangential kinetic energy is of the same order of 
magnitude as the energy in the transversal modes. The 
geometry of the wave guide plays no role at this level. 
On the other hand, in the literature concerned primarily 
with the effects of the geometry of constraint manifolds 
[iHS] on the effective Hamiltonian, it is assumed that the 
constraining forces are "constant" along the constraint 
manifold. This is because the geometric effects are much 
smaller and would be dominated by the adiabatic po- 
tential otherwise. It is thus assumed that the tangential 
kinetic energy is of the same small magnitude as the geo- 
metric effects and thus much smaller than the transversal 
energies. 

In this paper we show how these two regimes are re- 
lated and derive an effective Hamiltonian valid on all 
interesting energy scales. It contains contributions from 
the adiabatic potential, from a generalized Berry connec- 
tion and from the intrinsic and extrinsic geometry of the 
constraint manifold. The derivation is based on super- 
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adiabatic perturbation theory and a mathematically rig- 
orous treatment of the problem is given in [¥j . We present 
our results first on a general and abstract level. However, 
there are several concrete applications which have moti- 
vated us and the many predecessor works, most notably 
molecular dynamics and adiabatic quantum wave guides. 
In Section |IV| we apply our results to adiabatic quan- 
tum wave guides and, in particular, obtain new results 
about global geometric effects in quantum wave circuits, 
i.e. closed wave guides. 



A. Qualitative discussion of the results 

Although the mathematical structure of the linear 
Schrodinger equation 
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-A* + T/*=:ffV, '^It^o^^L^A^dr), (1) 



is quite simple, in many cases the high dimension of the 
underlying configuration space A makes even a numerical 
solution impossible. Therefore it is important to identify 
situations where the dimension can be reduced by ap- 
proximating the solutions of the original equation (fl]) on 
the high dimensional configuration space A by solutions 
of an effective equation 
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on a lower dimensional configuration space C. The fac- 
tor C*^ allows for the possibility of additional internal 
degrees of freedom in the effective description. 

A famous example for such a reduction is the time- 
dependent Born-Oppenheimer approximation: Due to 
the small ratio ^^^^^ of the mass md of an electron and the 
mass TTinu of a typical nucleus, the molecular Schrodinger 
equation 
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V^, (3) 



on the full configurations space IR'^("+™) = ^ of electrons 
and nuclei, may be approximated by an equation 

i^*V^ = - s^^-^ + ^ci^ , V'|t=o e l2(m3", dx) , 

on the lower dimensional configuration space M'^"^C of 
the nuclei only. In this case the interaction V{x, y) of 
all particles is replaced by an electronic energy surface 
-Eel (a;), which serves as an effective potential for the dy- 
namics of the nuclei. The assumption here is that the 
electrons remain in an eigenstate of the electronic Hamil- 
tonian He{x) — ~^/S.y + V{x,y) corresponding to the 
eigenvalue Ec\{x). This assumption is typically satisfied, 
since the light electrons move fast compared to the heavy 
nuclei and thus the electronic state adjusts adiabatically 
to the slow motion of the nuclei. This is an example 
of adiabatic decoupling where the reduction in the size 
of the effective configuration space stems from different 
masses in the system. 

A physically different but mathematically similar situa- 
tion where such a dimensional reduction is possible are 
constrained mechanical systems. In these systems strong 
forces effectively constrain the system to remain in the 
vicinity of a submanifold C of the configuration space A. 
For classical Hamiltonian systems on a Riemannian man- 
ifold {A, G) there is a straight forward mathematical re- 
duction procedure. One just restricts the Hamilton fimc- 
tion to C's cotangent bundle T*C by embedding T*C into 
T*A via the metric G and then studies the induced dy- 
namics on T*C. For quantum systems Dirac [5" proposed 
to quantize the restricted classical Hamiltonian system 
on the submanifold following an 'intrinsic' quantization 
procedure. However, for curved submanifolds C there is 
no unique quantization procedure. One natural guess 
would be an effective Hamiltonian Hcs in (pi) of the form 



ffn 



-Ac + V^lc, 



(4) 



where Ac is the Laplace-Beltrami operator on C with 
respect to the induced metric and V^|c is the restriction 
of the potential V : A^Rto C. 

To justify or invalidate the above procedures from first 
principles, one needs to model the constraining forces 
within the dynamics (II]) on the full space A. This is done 
by adding a localizing part to the potential V. Then one 
analyzes the behavior of solutions of (IT]) in the asymptotic 
limit where the constraining forces become very strong 
and tries to extract a limiting equation on C. This limit 
of strong confining forces has been studied in classical 
mechanics and in quantum mechanics many times in the 
literature. 

The classical case was first investigated by Rubin and 
Ungar [6;, who found that the effective Hamiltonian for 
the motion on the constrained manifold contains an extra 
potential that accounts for the energy contained in the 
normal oscillations. The quantum mechanical analouge 
of this extra potential is the adiabatic potential. The 
intrinsic geometry of the submanifold only appears in the 



definition of the kinetic energy ^g{p,p), its embedding 
into the ambient space A plays no role. 
On the other hand, for the quantum mechanical case 
Marcus [7] and later on Jensen and Koppe [H] and Da 
Costa [S] pointed out that the limiting quantum Hamil- 
tonian contains a potential term, the geometric potential, 
that depends on the embedding of the submanifold C 
into the ambient space A. But these statements (like 
the more refined results by Froese-Herbst [T], Maraner 
"2] and Mitchell [3]) require that the constraining poten- 
tial is the same at each point on the constraint manifold. 
The reason behind this assumption is that in the limit of 
strong confinement the adiabatic potential is much larger 
(by two orders in the adiabatic parameter) than the ge- 
ometric potential. For the geometric potential to be of 
leading order one must thus assume that the tangential 
kinetic energy is of the same small order. Then one ends 
up in the situation where the energy in the transversal 
modes is much larger than the typical tangential ener- 
gies and where, by assumption, any transfer of energy 
between transversal and tangential modes is suppressed. 
In conclusion, the effective Hamiltonian obtained in this 
way describes the constrained system only for very small 
energies and under very restrictive assumptions on the 
confining potential. Note that in many important appli- 
cations the assumption of a constant confining potential 
is violated. For example for the reaction paths of molecu- 
lar reactions, the valleys vary in shape depending on the 
configuration of the nuclei. 

In this work we present a general result concerning the 
precise form of the limiting dynamics |2]) on an arbitrary 
constraint manifold C starting from (JTF on the ambient 
space A with a strongly confining potential V. The most 
important new aspect of our result is that we allow for 
confining potentials that vary in shape and for solutions 
with normal and tangential energies of the same order 
and, at the same time, capture the effects of geometry. 
As a consequence, our effective Hamiltonian on the con- 
straint manifold has a richer structure than earlier re- 
sults and resembles, at leading order, the results from 
classical mechanics. However, similar to the hierarchic 
structure of the spectrum of molecules, with electronic, 
vibrational and rotational levels, now different geometric 
effects appear in the higher order corrections to the effec- 
tive Hamiltonian. We note that in the limit of small tan- 
gential energies and under the same restrictive assump- 
tions on the confining potential we recover the limiting 
dynamics by Mitchell [3]. 

The key observation for our analysis is that the problem 
is an adiabatic limit and has, at least locally, a struc- 
ture similar to the Born-Oppenheimer approximation in 
molecular dynamics. In particular, we transfer ideas from 
adiabatic perturbation theory, which were developed 
on a rigorous level by Nenciu-Martinez-Sordoni [TDHT^ 
and Panati-Spohn-Teufel in |13H15j and independently 
on a theoretical physics level by Belov-Dobrokhotov- 
Tudorovskiy in (THj , to a non-flat geometry. We note that 
the adiabatic nature of the problem was observed many 



times before in the physics Hterature, e.g. in the context 
of adiabatic quantum wave guides and thin films [171 HI] ■ 
But the only work considering constraint manifolds with 
general geometries in quantum mechanics from this point 
of view so far is }19j , where only the leading order dynam- 
ics of localized semiclassical wave packets is analyzed and 
effects of geometry or geometric phases play no role. We 
thus believe that our effective equations have not been 
derived before, neither on a mathematical nor on a the- 
oretical physics level. 



B. The scaling explained in a simple example 

Before we describe the general setup, it is instructive to 
first explain the scaling and the different energy scales 
within the simple example of a straight quantum wave 
guide in two dimensions. Let x be the coordinate in the 
direction of propagation and y the transversal direction. 
Saying that the potential V{x,y) is (at least partially) 
confining in the y direction just means, that the nor- 
mal or transverse Hamiltonian H,^(x) :— —A.y + V{x,y) 
has some eigenvalues Ej{x) with localized eigenfunctions 
ipj{x,y), the constrained normal modes. For a sketch 
of such a potential see Figure [TJa). Now we would 
like to implement the asymptotic limit of strong con- 
finement in such a way, that the eigenfunctions of the 
scaled Hamiltonian H^{x) become localized on a length 
scale of order e ^ 1. This is done by scaling the poten- 
tial V^{x,y) := V{x,y/e), which yields restoring forces 
of order e~^. However, localization on a scale of order e 
leads to kinetic energies of order e~^. So in order to see 
localization one has to increase not only the forces but 
also the potential energies to the same level by putting 
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'V{x,y/e). 



Then the normal energies and eigenfunctions are just 
E^{x) — e~'^Ej{x) and (p'j{x,y) = e~^ipj{x,y/e). The 
full Hamiltonian becomes 



H' = -A, - A, 



'V{x,y/e). 



In order to understand the asymptotic limit e — >■ it is 
more convenient to rescale units of energy in such a way 
that the transverse energies are of order one again, i.e. 
to look at 



H' 



e^H' 



-e'^.,~e'^y + V{x,y/e). (5) 



Changing units of length in the transverse direction to 
y = y/e finally leads to the form of the Hamiltonian 



H, := e^H^ = -e^A, - A,. 



V{x,y). 



(6) 



for which the normal eigenfunctions tpj are independent 
of e and the physical meaning of the asymptotic e — >■ 
is most apparent. The limit of strong confinement re- 
ally corresponds to the situation where the transversal 
modes are quantized with gaps of order one, while in the 



tangential direction the behavior is semiclassical and, in 
particular, the level spacing is of order e^. Here it is easy 
to guess the leading order effective Hamiltonian for the 
constrained system: on the subspace of wave functions of 
the form '^'^(x,y) — Lpj{x^y)'ip^{x) the Hamiltonian acts 
as 

i/e*"(a^,y) = {-e^/^r.-^y + V{x,y))^j{x,yW{x) 
= ^j{x,y) [(-e^A, -f ^j(x)) ^'{x)] 
-2{e\7^(Pj{x,y)){e\7^iP%x)) 
~e^A,^j{x,y)^'{x). (7) 

Defining the effective Hamiltonian by projecting back 
onto this subspace via P{x) := \ipj{x,-)){(pj{x,-)\ and 
integrating out y, one finds 

H!b = HsV.,-is{^,ix)\V,^,{x))f+E,{x) (8) 

+ e^{{V^ip,{x)\V,ip,ix)) - \{ip,{x)\V^ipjix))\^). 

Here we see how the transversal eigenvalue Ej{x) en- 
ters as an effective potential, the adiabatic potential, at 
leading order. E.g., the constraining potential sketched 
in Figure ffja) leads to an attractive effective potential 
sketched in Figure [2J The two energy scales referred to 
in the previous section now correspond to the following 
situations: if one assumes "small" tangential energies, 
i.e. {ip^\ - e'^A^tp") = 0{e^), then all terms in Q but 
the term involving the adiabatic potential Ej(x) are of 
order e^. Thus the latter must be either constant or the 
kinetic energies will also become 0(1) under the time evo- 
lution. Since it turns out that the geometric potential in 
the case of non-straight wave guides is also of order £^, 
this explains why all authors interested in geometric ef- 
fects up to now assumed Ej{x) = const. 
However, the natural scaling is to allow for tangen- 
tial states ip"^ with kinetic energies of order one, i.e. 
(-0^1 — s'^A^ip^) = 0(1). Then all energies in the sys- 
tem are of the same order and exchange of normal and 
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FIG. 1: In (a) we plotted a potential for a waveguide which 
widens near x — 0. The widening lowers the energy of normal 
modes and thus produces an attractive effective potential for 
the motion in a;-direction. In (b) the modulus of the ground 
state wave function is sketched. Its variation in a;-direction 
is slower than in y-direction, but its tangential derivatives al- 
ready grow in e. In (c) the modulus of an excited state with 
energy of order one above the ground state is sketched. Its 
variation in a;-direction is on the same scale as the confine- 
ment, i.e. it oscillates on a scale of order e. Thus any analysis 
assuming bounded tangential derivatives of the solutions will 
be restricted to confining potentials with constant profile. 




Eo{x) 



0(1) 



highly nontrivial. Roughly speaking, for times of order 
one the problem is just standard time-dependent pertur- 
bation theory. For times of order e^^ one can use the 
ideas underlying the standard proof of the adiabatic the- 
orem of quantum mechanics, see Subsection III A For 
longer times, however, one has to use "super" -adiabatic, 
i.e. higher order adiabatic perturbation theory, see Sub- 
section unci 



II. THE ADIABATIC STRUCTURE 



FIG. 2: The curves Eo{x) and Ei{x) are sketches of the low- 
est normal eigenvalues for a waveguide potential as depicted 
in Figure flja) . On the vertical axis the spectrum of H^ is 
drawn: one expects spectral bands starting at the minima eo 
and ei of the effective potentials with level spacing of order e^. 
The continuum edge E is determined by the threshold of Eo . 
Eigenstates in the lower shaded region vary on a y^-scale in 
the a;-direction as indicated in Figurellfb). Eigenstates in the 
upper shaded region with energies of order one above eo have 
e-oscillations in the a;-direction as indicated in Figure [T][c). 



tangential energies may occur. In particular, the tangen- 
tial momentum operator —ie'Vx must be treated as being 
of order one despite the factor e. This is the situation we 
will consider in the following. 

In the Figures [T] and [2] we sketch the situation for a 
simple waveguide in a region where it widens slightly. 
Wave functions with tangential energies of order e like 
in Figure [l[b) yield the low lying part of the spectrum. 
General states with finite energy above the ground state, 
which include all states propagating through the wave 
guide, have tangential energies of order one and thus e- 
oscillations in a;-direction, as indicated in Figure fl|c). 
When the confining potential depends on x, there are, in 
general, no solutions with tangential kinetic energies of 
order e^ . We also mention an extensive discussion of en- 
ergy scales from a slightly different point of view in |17j . 
Before we explain the general model, it is instructive to 
mention two important points on the level of this simple 
model. First of all one might want to add an "external 
potential" W{x, y) which does not contribute to the con- 
finement and thus is not scaled. We will allow for such an 
external potential with W{x, y) = 0(1). Note that in the 
previous works focussing on geometry [Ij [3] it was added 
on the small energy scale, i.e. W{x,y) = 0{e^). The 
second remark is that with the energy scales there come 
also time scales. The time scale on which solutions with 
(-0^1 — e^^x^'^) — C(l) propagate distances of order one 
are times of order e~^ . This is because kinetic energies 
of order one for particles with "mass" of order £~^ yield 
velocities of order e. The small energy solutions with 
(V'^l — e'^Axip'^) = 0{e'^) propagate even slower, so here 
the natural time scale are times of order e~^. The best 
results we can prove hold for even longer times, namely 
for times almost up to order e~^. Controlling the adia- 
batic decoupling for such long times makes the problem 



Here we first discuss in detail the model we consider and 
the assumptions involved. In the second subsection we 
introduce a horizontal momentum operator, the geomet- 
ric generalization of —is'Vx in the previous section, which 
will play a crucial role in our results. Then we reveal the 
formal similarity with the Born-Oppenheimer approxi- 
mation, before we explain the resulting adiabatic struc- 
ture of the problem in the last subsection. 



A. Description of the model 

Let (A, G) be a Ricmannian manifold of dimension d + k 
and C (Z Aa. smooth submanifold of dimension d without 
boundary and equipped with the induced metric g = G\c- 
We consider the Schrodinger equation on A with a po- 
tential V^ : A ^ M. that localizes all states from a certain 
subspace of LP' {A) close to C for small e, which will be 
made precise below. 

So we want to start with fixed manifolds A and C and 
assume that the constraining potential V^ grows fast in 
the directions normal to C (strong restoring forces) while 
the variation along C is of order one (bounded tangential 
forces). This means we want to assume that 

• normal derivatives of V^ are of order £~^ , 

• tangential derivatives of V^ are of order 1, 

• all derivatives of the metric G are of order 1. 



As explained in Subsection I B localization in the normal 
direction on a scale of order e produces oscillations of 
order e~^ in the tangential directions, too. 
When we introduce local coordinates x = (xi)i=i...._d in a 
neighborhood of g e C and coordinates N — {Na)a=i,...,k 
for the normal directions, the assumptions made above 
correspond, by the same reasoning that leads to ([5]) in 
Subsection |IB[ to the Schrodinger equation 



ia*^ 



- e^Ag*^ + Vc{x, iV/e)*^ + W{y, N)'^' (9) 



where A^ is the Laplace-Beltrami operator associated 
to {A, G) and the non-constraining potential W may de- 
scribe external forces. Here the upper index e at VP means 
that we look at solutions with oscillations of order e~^ 
and the e^ in front of the kinetic energy ensures that these 
solutions have kinetic energies of order 1. For small e at 



least some solutions of this equation concentrate close to 
the submanifold C. Therefore one expects that an effec- 
tive Schrodinger equation on C may be derived such that 
solutions ?Ae {t) of the effective equation approximate the 
solutions ^"^ (i) of the full equation in a suitable way. 

The scaling of the potential described in (|9| depends on 
the choice of coordinates and cannot be implemented 
globally so naively. It just serves as a motivation for 
the following. In order to be able to implement a sim- 
ilar scaling globally we assume that the submanifold C 
has a tubular neighbourhood B of fixed diameter (5 > 0. 
Within B it makes sense to speak of large derivatives of 
V^ with respect to the distance to C. More precisely, B 
can now be mapped to the (5-neighbourhood Bs of the 
zero section in the normal bundle NC. On NC the scal- 
ing of the potential as in (l9| can be realized due to its 
linear structure. Moreover, for e much smaller than 5 
all solutions below an arbitrary finite energy lie in Bs/2 
up to errors bounded by any power of e. Therefore it 
is possible to work completely on the normal bundle by 
constructing a diffeomorphism <I> : NC — >■ B and choosing 
a metric g on NC such that $ is an isometry on Bs/2- 

To avoid all regularity problems we make the following 
assumption. 

Assumption 1: The injectivity radii of A and C are 
strictly positive and all curvatures as well as their deriva- 
tives of arbitrary order are globally bounded. Further- 
more, V : NC -^W is smooth and bounded and arbitrary 
derivatives of V are also globally bounded. 

In particular, this implies that C and NC may be covered 
by coordinate neighborhoods such that for some K G N 
not more than K of them overlap at each point. This 
allows us to do all estimates in local coordinates. 

Our goal is now to find approximate solutions of the 
Schrodinger equation 

i^t^^ = -£2^g*^ + Vciq, N/e)"^' + W{y, N)^'^ 

on H = L'^{NC,d]l), where g is the puUback of G via 
the diffeomorphism $ on Bs/2j suitably extended outside, 
and dJI denotes the measure associated to g. As expalined 
in the introduction it is helpful to rescale the normal 
coordinates to n = N/e. Then the equation reads 

i9t*^ = -e^{AgY^^+Vciq,n)'^^+W{y,en)^^ =: H,^^ , 

where (A^)^ is the accordingly rescaled Laplacian, whose 
expansion in e is calculated in the appendices A2-A4. 



B. The horizontal connection and the 
corresponding Laplacian 

Since we will think of the functions on NC as map- 
pings from C to the functions on the fibers, the follow- 
ing objects will play a crucial role. Consider the bundle 



£f := {(g, ip)\q G C, ifi G L'^{NqC)} over C which is ob- 
tained when the fibers NqC of the normal bundle are 
replaced with L?[NqC) and the bundle structure of NC 
is lifted by lifting the action of SO(fc) on the fibers to 
rotation of functions. We denote the set of all smooth 
sections of a hermitian bundle £ by r(£). 
For (f G r(£f ) the horizontal connection V'' is defined by 



d_ 
ds 



s=0 



(p{w{s),v{s)), 



(10) 



(V^^)(g,n) 

where r € TgC and (w, v) e C^hl, l],NC) with 

w{0) = q, w(0) = T, & u(0) = n, V^v = 0. (11) 

Furthermore, Ah is the bundle Laplacian associated 
to V'^, i.e. defined by 

/ / ip*Ai,^pdndn ^ - f f g'^^f^'^];ipdndfi 

JC JNgC JC J NgC 

where g'^^ is the inverse of the metric tensor gij. Here 
and in the sequel we use the abstract index formalism 
including the convention that one sums over repeated 
indices. Moreover, we will consistently use latin indices 
«,j, .. running from 1 to d for coordinates on C, greek 
indices a, /3, . . . running from d-|-l to d+Zc for the normal 
coordinates, and latin indices a,b,.. running from 1 to 
d -\- k for coordinates on the full normal bundle. 

To obtain local expressions for these objects we fix g e C 
and choose geodesic coordinate fields {dxi}i=i,....d on an 
open neighborhood f2 of g and an orthonormal trivializing 
frame {va}a=i,...,k of Nft. We define the connection co- 
efficients u!^^ of the normal connection by V^i^Q, = u^^i^j. 
Then the horizontal connection is given by 



V^(y5((7,n) = d^ip{q,n)~uj]^n°'djip{q,n), (12) 

as was already shown by Mitchell [3j , and it holds 

Ah(^ = fi-^d, - ujfpn^ds)pig''{dj - ojj^n'^d^)^ (13) 

with fi := det gij. The latter directly follows from the 
former and the definition of A^. To obtain the former 
equation we note that for a normal vector field v — n'^Va 
over C it holds 



{VivV = d.n-^+uj-l^n^. 



(14) 



Now let (w,w) G C^[[~l,l],NC) be as in ^. Then by 
definition of V^ we have 

= diip{q,n) + {din'')d^(p{q,n) 
= dt(p{q,n) - uj^^n"'d^ip{q,n) 



where we used ( 14 ) and the choice of the curve v in the 
last step. 



C. The splitting of the Laplace-Beltrami operator 

The basic idea for deriving an effective equation on the 
submanifold C is to spht the Laplace-Beltrami operator 
on NC at leading order into a horizontal and a normal 
part relative to C, similar to the splitting — A^ — A^ in 



the simple example of Section IB To make this precise, 
first note that by construction at any point on the zero 
section of NC (which we identify with C in the following) 
the tangent space splits into two orthogonal subspaces, 
one tangent to C and one tangent to the fibre. Hence the 
metric tensor g and with it also the Laplace-Beltrami 
operator on {NC,g) splits into a sum 



A^ 



^N, 



where Ag is the Laplace-Beltrami operator on C and An 
is the euclidean Laplacian in the fibers NqC = M.'' of the 
normal bundle. We note that A„ — Ah on functions that 



are constant on the fibers by ( 13 ). We will show that also 



away from C, i.e. globally on NC, we can approximately 
split Ag into a horizontal part, given by Ah, and the 
Laplacian in the fibre Ajv- The error grows linear with 
the distance |iV| to C Then the rescaling of the normal 
coordinates to n = N/e yields that 

H, = ~e^Ai,-A^ + Vc{q,n) + Wiq,en)+0{e\n\). (15) 

This operator has the same form as the Hamiltonian pi), 
which is the starting point for the time-dependent Born- 
Oppenheimer approximation, or the operator ([6| of our 
simple wave guide example. This suggests that also in the 
general situation considered here adiabatic decoupling is 
the mechanism that yields effective Hamiltonians on C. 
We now explain in more detail how to achieve the above 
splitting of the Laplacian. An important step is to turn 
the measure on NC into product form. To do so we define 



Mp : L^{NC,dJI) 



L^{NC,dNdii), 



* H^ Mp^^ := p 



*. 



where dN denotes Lebesgue measure on the fibers NqC = 
R*^ and p = ^ J^^j^ is the density of the original measure 
with respect to the product measure on NC. It is well- 
known that the unitary transformation of our Hamilto- 
nian with Mp leads to the occurence of a purely geometric 
extra potential 

Vp = - \r'{da In p) db In p + i Ag In p. 

More precisely, it holds that 

M;{~Ag)Mp = -Ah - Aw + Vp{q,N) + 0{\N\) , 

(16) 
which is shown in the second to fourth appendix. There- 
fore after application of the unitary transformation Mp 
and a Taylor expansion of W the rescaled Hamiltonian 
H^ is of the following form close to C: 



H, 



-e^Ah - A,, 



Vciq.n) + Wiq,0) 
^e^Vp{q,sn) + 0{e\n\) 



We note that — e^Ah is of order 1 on functions with os- 
cillations of order e~^. So the extra potential does not 
play a role for the leading order of the horizontal dynam- 
ics, unless the tangential kinetic energies are assumed to 
be small. Finally, it should be kept in mind that the 
remaining error term is small only when it is applied to 
functions that decay fast in the normal directions. 



D. Adiabatic decoupling 

Next we explain the principle of adiabatic decoupling in 
detail. For any q G C we define the fiber Hamiltonian 

Ht{q) = -A„ + V,{q,n) + W{q,0) 

on the Sobolev space W'^'^{NqC,dn) C L^{NqC,dn). We 
consider a g-dependent family of eigenvalues Ef{q) of 
multiplicity M, called an energy band in the sequel, 
and an associated family of normalized eigenfunctions 



,M' 



H,{q)^f{q,-) = E,[q)vi{q,-)- 

By definition of Hf it holds i/^ — H{ — e^A^ - 
on states that decay fast enough. Then states in 



(17) 
0{e) 



Vo ■■= Wi'{q,n)^Pjiq)\ijj e L^iCg)} C L^iNC) 

are approximately invariant under the dynamics e~'* = 
for times of order 1. This is due to the fact that the asso- 
ciated projector Pq defined by Po{q) '■— {\(Pf){iPf\){q) 
is a spectral projection of H{ and so we know that 
[H{, Po] = 0, [^f , Pa] = 0, and HfPa = EfP^. Hence, 



[H„Po] - [-£2Ah,Po] + Oie) = 0{e). 



(18) 



More precisely, the solution of the full Schrodinger equa- 
tion with initial value ^"^jt^o = "^j v/ satisfies that 

^'{t,q) = ip,\q,nWjit,q) + Oie\t\) , 

where 'ipj{t,q) solves the following effective Schrodinger 
equation on C: 



idtTp'y{t,q) 



-e^Aglb'j{t,q) + Et{q)rj{t,q). (19) 



It is well-known that an equation of the form ( 19 1 does 



show interesting behavior only on the semiclassical time 
scale s = t/e. The adiabatic principle, however, suggests 
that Vo may be expected to be invariant for such and even 
much longer times, if the energy band (-E'f)g is separated 
by a gap from the rest of the spectrum. Therefore we 
also assume the following. 

Assumption 2: For all q E C the fiber Hamiltonian 
Hi{q) has an eigenvalue E{{q) of multiplicity M such that 

inf dist(^f(g), spec{Hi{q) \Ei{q)))> Cgap > 0. 



In addition, there is a family of normalized eigenfunc- 
tions {ff{Q))j^i jyj which is globally smooth in q (in 
particular, the corresponding eigenspace bundle is trivi- 
alizable) and satisfies 



A. Effective dynamics for times of order e ^ 

The analogue of the adiabatic theorem, which provides 
effective dynamics for times of order £~^, reads as follows: 



supgec lle^''^"^/!! <oo, 



for (n) := -^/l + |np, some Ap > 0, and all J. 

An assumption about the decay is necessary because the 
error in the splitting — Ag — — e^A^ — A„ + 0{e) is only 
small when applied to functions that decay fast enough, 
as was explained above. However, in lots of cases the 
decay is implied by the gap condition, in particular, for 
Ef below the continuous spectrum of Hf. The assumption 
about triviality is necessary to get an effective equation 
on -L^(C, d^) ® C*-'^. If we dropped it, we would end up 
with an equation on a non-trivial rank-M bundle over 
L^(C,d/x), which would complicate things quite a bit. 



Theorem 1 Fix E-n-^a^ < oo and denote by x the charac- 
teristic function of (— oo,£'max]- Lst the energy band E{ 
and the family of normalized eigenf unctions (<p'f )j=i,...,m 
be as in Assumption^ 
Define the operator Uq-.To^ i^(C, d/i) (g> £^' by 

and H^^ :— UqPqHi^PqUq . Then there is a C < oo such 
that for all e small enough 

(e-'^^* - t/*e-'<'*[/o) Pox{He)\\ < Ce{l + e\t\). 

(20) 
Up to terms of order e^ the first-order effective Hamilto- 
nian H^g is given by 



III. MAIN RESULTS 

Having revealed the adiabatic structure of the constrain- 
ing Hamiltonian H^ in the preceding section we have two 
different techniques at hand in order to deduce results 
about effective dynamics. 

On the one hand, it is possible to derive an analogue of 
the standard adiabatic theorem of quantum mechanics 
in order to show that the subspace Vq is invariant under 
-ffj for times of order e~^ up to errors of order e. This 
is analogous to the approach used in |20| in the context 
of the Born-Oppenheimer approximation and will be car- 
ried out in the first part of this section. It leads to the 
occurrence of a Berry connection that will be investigated 
in the second subsection. 

In order to get a better approximation of the spectrum 
and/or to go to longer time scales for the dynamics, the 
usual adiabatic technique relying on cancellation of er- 
rors due to oscillations is no longer practicable. How- 
ever, the general machinery of adiabatic perturbation 
theory, developed by Nenciu-Martinez-Sordoni in pDHT^ 
and Panati-Spohn-Teufel in [T31II1] and reviewed in [T5] . 
allows to construct super-adiabatic subspaces that are 
invariant for times of order e^"i up to errors of order 
e"^ for arbitrary ni,n2 € N. More precisely, it allows 
to construct a projector P^ which projects to a subspace 
Ve close to Vo and satisfies [H^,P^] = ©(e™) for any 
TO > 1. Adiabatic perturbation theory was adapted to 
constrained quantum systems in [5j. For technical rea- 
sons it could be made rigorous only for to, < 3. The case 
m — 3 seems enough for all applications though. Before 
we discuss the resulting effective Hamiltonian and the ap- 
proximation of bound states in the last two subsections, 
we explain the construction of P^ in the third subsection. 



(^^|i?Wv^)c = 



{g%"p'I,rKPiijrL + vj/rj^j) d/. 



with 



PcSj 



ie<5''^a, -e(^/|iVX> 






ffeff 



IJ 



e2IlL^"(^f^|n>/) 



V/„' = Ef6'-' + eid^W)n=o{^iK^f) 



where II is the second fundamental form (see Appendix 1 
for the definition) , ( ■ | • )c is the scalar product on 
L?{C,d^), and {■]■) is the scalar product on L'^{M.^ , dN) . 

Via the operator Uq it is, hence, possible to obtain ap- 
proximate solutions of the original equation from the 
solutions of the effective equation. We point out that 
PoXiHe) both cuts off high energies and produces initial 
states in Vq. However, the cutoff energy -Emax is arbitrary 
and, in particular, independent of e. It is only needed in 
order to get a uniform error bound, since for larger tan- 
gential energies the adiabatic decoupling becomes worse. 
Physically this is expected, since large tangential energies 
correspond to large tangential velocities and the sepa- 
ration of time-scales for the normal and the tangential 
motion, which adiabatic decoupling is based on, breaks 
down for large velocities. 

The effective Hamiltonian may be calculated using stan- 
dard perturbation theory which is done in the last ap- 
pendix. However, to verify that it yields effective dynam- 
ics on the relevant time scale t = s/e, i.e. to prove (20 1, 
an additional adiabatic argument is needed. To make this 
clear we notice that the usual perturbative argument only 
yields an error of order 1 for times of order e~^: Using 



that UqUq ~ Pq and UqUq ~ 1 we have 
fe-'^"* - C/*e-'^''^°^^^°^o*C/n') Pn 



(21) 



-iH^t 



iH^t 



d^ iHcSTj* -iUoH^UgS 
ds 



o^Uods 



e'"^n[H„Po]USc 



-lUaH^UXs 



o'Uods, 



which is of order e|i| by (18) but cannot directly be seen 
to be smaU for times of order e~^ . However, adapting the 
calculation in the derivation of the standard adiabatic 
theorem (see e.g. [15]) shows that the integrand is, up to 
errors of order e^, the time derivative of 

e'^-^ {RHAHe,Po] ~ [Po,H,]Rh,) C/o*e-'^°^-^o^C/o, 

where Rhi ~ P(t i-^f ^ -^f)~~'^ -^o' ^^ ^^^ reduced resolvent. 
Therefore the time integral of this term yields an error 
of order e independent of t, which shows that the whole 
error is, indeed, only of order e{l + e\t\). 

For times of order e~^ the corrections of order e yield 
relevant contributions: 

• The corrected momentum operator p'^^ is a Berry 
connection on the C -bundle over C where the ef- 
fective wave function takes its values and there- 
fore may give rise to topological and/or geometric 
phases (see the next subsection). 

• In general, all the corrections couple the effective 
internal degrees of freedom. If they, however, mutu- 
ally commute, simultaneous diagonalization allows 
to split the effective C^^-bundle into M effective 
bundles of rank 1 locally. 

• If the center of mass of (pj(pf{q) lies on C for all 
q, then {(pf \ n°'ipf) = 0. So both the correction to 
g^g and to Ves vanish in this case. In particular, a 
local splitting into M effective bundles of rank 1 is 
possible in this case. 



B. The curvature of the Berry connection 



In this section we take a closer look at the induced Berry 
connection p^^ = —iedxS^'^ — e{ipf\iW^LPf) that occurs in 
the effective Hamiltonian (see Theorem [I]). 

For M = 1, i.e. if the energy Ef is non-degenerate, it is 
simply a [/(l)-connection that effects the dynamics sim- 
ilar to the vector potential of a magnetic field. If its 
curvature (the analogue of the magnetic field) is zero, 
one can achieve, at least locally, poff' = —^sdx by choos- 
ing a proper gauge, i.e. by choosing proper eigenfunc- 
tions (pf{x). But such a gauge might not exist globally 
and effects analogous to the Aharonov-Bohm effect may 
occur. In Section |IVB| we give an example for a closed 
quantum wave guide without any external magnetic fields 



in which such an effect occurs purely due to the geometry 
of the wave guide. 

If the curvature of Pcff is non-zero, it does even locally 
change the dynamics at order e. In [21] Maraner identi- 
fied the curvature as the origin of roto-vibrational cou- 
plings in simple molecular models. Moreover, further im- 
portant effects are to be expected which are known for 
Berry connections from different areas: On the one hand, 
the anomalous velocity term in the semiclassical model 
for electrons in crystalline solids also stems from the cur- 
vature of a Berry connection, see [55J |53| . On the other 
hand, in the Born-Oppenheimer approximation the Berry 
connection term in the effective dynamics exactly cancels 
the effect of an external magnetic field on the nuclei, see 
e.g. [20 . Neglecting the Berry term would lead to wrong 
physics: in the Born-Oppenheimer approximation a neu- 
tral molecule would suddenly react to the Lorentz force. 

In the rest of this section we show how to obtain the fol- 
lowing formula for the curvature of the Berry connection. 
To state it, we &x q € C and choose again normal coordi- 
nate fields {di}i^i_,,,^d on an open neighborhood fi of q. 
Then it holds [d„dj]^0. 

Proposition 1 V°^ :— ipoff is a metric connection on 
the rank-M bundle over C where the effective wave func- 
tion takes its values. Its curvature vanishes for M = 1 
and otherwise is given by 






cffyvcff 



VfV 



cffyycff 



Vf V 






= (vfvf-vfvf)^^ 

= -e'{^UK^]n'^d,^^){q) 



((Vh^f^|V>/)-(V>f^|Vl'(p/))(g) 



(^f^lv^VfX^flv'v/) 



{^i\ 



Vj(^f 



K\ 



{v\ 



K 



^Wi) 



(9) 



where R^ is the curvature of the normal connection (de- 
fined in the first appendix) and ( • | ■ ) is the scalar product 
onL^{R'',dN). 

An analogue expression was derived by Mitchell in [3] 
in the special case where (pf is independent of q up to 
twisting. It was not realized that it always vanishes for 
M = 1 though. 

We will need that the connection V*^, which the nor- 
mal connection induces on the bundle of functions 
over the normal fibers, is metric, i.e. dj{ipi\ip2)'Hf — 
{W^(pi\ip2)'Hi + {'Pi\'^]j'P2)'Hn ^^'^ t^^t its curvature is 
given by 



R!^ :- Vj^V^ 



V,^V^ 



[d.M,] 



R«]""a,. 



(22) 



Since the normal connection is metric, its connection co- 
efficients Luf^ are anti-symmetric in a and f3. So integra- 
tion by parts yields 

{uj^^n"djipi\ip2){q) + {(pi\uj]^n°'djip2)iq) = 0. 



Therefore we have 

dj{ipi\^2) = {dj^i\ip2) + {ifi\djif2) 

+ ((^i|(9j-a;7„n°5^)^2) 

To compute the curvature of V^ we notice that a simple 
calculation yields 

Using the commutator identity 

we obtain that 



Then (22) yields 



R''. 



i9x^^]a - 9x^,0.+ ^t^]p - u:^j^u:Ji^)n°'dn-, 






which was the claim. With this we can compute the cur- 
vature of the effective Berry connection. It is not difhcult 
to verify that V^* is indeed a connection. Since V*^ is 
metric, we have that 



Thus the correction in V°^ is anti-hermitean. Hence, for 

all ^1,^-2 :C^C*^ 

ea-,(i/'i • -02) = (ec'jV'i) • 0^2 + V'l • {£djip2) 

which means that V^ is metric. Furthermore, this en- 
tails that the correction in V°^ is purely imaginary for 
M = 1. Since ipf can be chosen real- valued for every 
q £ C, which follows from H{ being real, we may gauge 
away the correction in an open neighborhood of any q. 
This implies that the curvature vanishes for M = 1. 
To compute the curvature of V"**" for M > 1 we calculate 



R/Jij 



£2 (V^^V^* 



3 « JlJ 



= e^d,{^l\V^^f) - a,(^f|V^(p/)) 



(^f|v^^f)(v.f |VW 



Using that V^ is metric we obtain 

-(<^f|V>f)(^f|VjV/) 



R7;; = -e^{^',\Rt;in-d,^i) 

-f£2((VM|V>/)-(V>^|Vl'V'/> 

+e2((^f^|v5Vf)(</.f|V>/) 



-(^f|V^(^f)(^f|V^V'/; 



which was to be shown. 



C. Construction of the superadiabatic subspace 

There are several motivations and ways for further im- 
proving the result formulated in Theorem [l] First of all 
one can aim at a better approximation, i.e. smaller error 
estimates. Next one can try to cover even longer time 
scales, i.e. times of order £~^ and beyond. These long 
time scales become relevant, e.g., when considering the 
propagation of states with tangential energies of order e^ 
in wave guides where the energy band E[ is constant on 
all of C, i.e. in the situation considered in earlier papers on 
geometric effects on constrained systems [IHSl [3T] . Last 
but not least one expects also the eigenvalues of the effec- 
tive Hamiltonian to be close to those of the full Hamilto- 
nian and that one can recover, at least in a certain energy 
range, all eigenvalues of the full Hamiltonian in this way. 
In order to achieve all three additional goals we show 
how to construct an effective Hamiltonian that is unitar- 
ily equivalent to the full Hamiltonian on a certain sub- 
space of the full Hilbert space up to errors of order e^. To 
this end we use adiabatic perturbation theory [13|. The 
strategy is to first associate a so called super-adiabatic 
subspace V^ with any energy band E{ satisfying Assump- 
tion [2J The associated projector P^ turns out to be 
uniquely fixed (up to terms of order e^) by the require- 
ment that it projects on Vq to leading order and that the 
commutator [H^,P^] is of order 0{£^). In a second step 
we construct a unitary f/e mapping the range of Pe to the 
Hilbert space of the constrained system L'^{C, dfj.) (8) C*'^. 
Then on the super-adiabatic subspace H^\-p^ = PsH^Pe 
up to terms of order e^. The effective Hamiltonian on 



(2) _ 
eff 



u,pm,pm: 



L^{C,dii) (g) C^ is now given by iJ, 

and solves all three problems mentioned above. 

We now explain this construction in detail. For the super- 
adiabatic projection we search for a bounded operator P^ 
with 

i) P,P, = P„ 

ii) P,-Po = 0{e), 

iii) [H,,P,]x{H,) - 0(e3). 

Property i) simply means that Pg is an orthogonal pro- 
jection, property ii) is the requirement to be close to the 
adiabatic projection Pq and iii) says that P^x(He)'H is 



invariant under the Hamiltonian H^ up to errors of or- 
der e^. 



Since we saw in ( 18 1 that [H^, Pq] — 0{e), it is consistent 
to make the ansatz for P^ as an expansion in e starting 
with Pq: 



P. 



Pq + ePi 



e'P2 



0{e^ 



We first construct P^ in a formal way ignoring prob- 
lems of boundedness. Afterwards we will explain how to 
obtain a well-defined projector and the associated uni- 
tary U^. We make the ansatz Pi := T^Pq + PqTi with 
Ti : H ^- "H to be determined. Using the expansion of 
Hg = Hq + eHi + 0{e^) from the fourth appendix and 
assuming that [Pi, — e^Ah -f -Ef] = 0(e) we have 



[H,, P,]/e = [Ho/e + i/i, Pq + ePi] + 0{e) 

= [Ho/e + Hi,Po] + [Ho,Pi]+Oie) 

= [-eAh + i/i, Po] + [Hf - Eu Pi] + 0{e) 

= (-eAh + i7i)Po - Po(-£Ah + Hi) 

+ {Hf - Ei)T*Po - PoT.iHi - El) + 0{s) . 

We have to choose Ti such that the first term is can- 
celled. Observing that the right hand side is off-diagonal 
with respect to Pq, we may multiply with Pq from the 
right and P^ := 1 — Pq from the left and vice versa to 
determine Pi. This leads to 

-(i/f-Pf)"'Po^([-£Ah,Po]+i/i)Po = Po^TrPo (23) 

and 

^Po{[Po,-eAi,]+Hi)P^{Ht~Efy^ = PqPiPo^, (24) 

where we have used th at t he op erat or Hf — Ef is invertible 
on Po^Hf. In view of (psl and ([24|, we define Ti by 



Ti := -Po([Po,-£Ah]+iJi)Pff.. 

+Pff,(heAh,Po]+i?i)Po (25) 

with Rh{ = P^[H{ — Sf) Pj-. Ti is anti-symmetric so 
that P(i) := Pq + ePi ^ Pq + e{T;PQ + PqTi) automat- 
ically satisfies condition i) for P^ up to first order: Due 
to P2 = Po 

p(i)p(i) ^ p(i)+ePo(Ti*+Ti)Po + C'(e2) 
= p(i)+0(£2). 

Moreover, it turns out that Pi satisfies the assumption 
[Pi, — e^Aji + E[] = 0{e) made above, too. 
In order to derive the form of the second order correction, 
we make the ansatz P2 = T^PoTi + T^Pq + P0P2 with 
some T2 : H — >■ "H. The anti-symmetric part of T2 is 
determined analogously with Ti just by calculating the 



commutator [P^,H^] up to second order and inverting 
H{ — E{. One ends up with 

(r2-T2*)/2 = -Po([P(l),ff(2)]/^.2)p^^. 

+RH,{[H('\P^'^]/e')Po 

with ij(2) := Ho + eHi + e^H2. Note that 
[ij(2)^p(i)]/£2 ^ 0(^) j^g i-Q ^Yie construction of P^^). 
The symmetric part is again determined by the first con- 
dition for Pe. Setting P^^' := P'^' + e'^P2 we have 

p(2)p(2) ^ p(2) ^ £2p^^p^p* ^ p* ^ j.^-^p^ _^ 0^^3)^ 

which forces r2* +T2 — —Pi Pi* in order to satisfy condi- 
tion i) up to second order. 

We note that Ti is quadratic in the momentum (and T2 
even quartic) and will therefore not be bounded on the 
full Hilbert space and thus neither P^. This is related 
to the well-known fact that for a quadratic dispersion 
relation adiabatic decoupling breaks down for momenta 
tending to infinity. The problem can be circumvented 
by cutting off high energies in the right place, which was 
carried out by Sordoni for the Born-Oppenheimer setting 
in [TT] and by Tenuta and Teufel for a model of non- 
relativistic QED in [21] . 

To do so we fix Pmax < 00. Since i/^ is bounded 
from below, P_ := mi(j{H^) is finite. We choose 
X e C(f(M, [0,1]) with x|(£;_-i,B+i] = 1 and suppx C 
(P_ - 2, P + 2]. Then we define 

P* := Po + (P(2) - Po)x{He) (26) 

+x(i?e)(P(')-Po)(l-x(i^s)) 

with x{Hg) defined via the spectral theorem. We empha- 
size that P/ is symmetric. 

It holds that P/ — Pq = 0(e) in the sense of bounded 
operators. That is why for e small enough a projector is 
obtained via the formula 



P. := 



27r 



iP} 



dz, 



(27) 



where r = {zGC||z — 1| = 1/2} is the positively ori- 
ented circle around 1 (see e.g. [25j). Denoting the as- 
sociated subspace by Ve we define a unitary mapping 
J/e : Pe ^ Pq by the so-called Sz.-Nagy formula: 



U, 



(PoPe + (l-Po)(l-Ps))(l-(Ps-Po)') '■ (28) 



Then Ug := UqUc yields an isometry between Ve and 
L^{C, d/i) (E) C*^. In 4] it is shown that P^, indeed, satis- 
fies i) to iii): 

Proposition 2 Fix Pmax < 00. For all e small enough 
Pj is an orthogonal projection and U^ is unitary. There 
are constants Ci such that 

• ||Pe-Po||£(«) < Cie, 

. \\[H,,P,]x{H,)\\ci'H.v{H^)) < C^e^ 

• \\{7iYP,{nY\\cm < C3 Vi,/GNo, (29) 
with X the characteristic function of (—00, E]. 
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The last estimate guarantees that the range of P^ consists 
of states decaying faster than any polynomial, which is 
necessary to use the expansion of i/^ obtained in the 
fourth appendix. 



D. 



Effective dynamics for times of order e 



By combining the results of the previous section with 
standard perturbation theory we can conclude that for 



i/^(2) ._ u^P.H^PeU; we have 



iH,t/ 



Ule-'^'^^'UA P.xiHe) 



<Ce''\t\ 



In the super-adiabatic setting no further adiabatic aver- 



aging is needed. This clearly improves (20) in the two 



ways anticipated: we get a better approximation and 
longer times. To get a simpler expression we can ap- 
proximate Pg and Ue by Po and C/q and find 



-iH^t 



U^e 



iH't t, 



«*C/o Pox(^e) <Ce{l + e'\t\), 



i.e. still a good approximation for long times on the adi- 

(2) 

abatic subspace Vo- However, we can not replace i?eff 
by H^g ^ UoPoHePqU^ without loosing a factor e in 
front of \t\ in the error. This is because the order e^ 
terms in the effective Hamiltonian are relevant for times 
of order £^^ and the expansion of the "naive" adiabatic 
Hamiltonian H^g yields incorrect second order terms. 

So we still have to provide the correct second order ex- 

(2) 

pansion of the effective Hamiltonian -ff^ff . Since the 
expression becomes quite complex and since we do not 
want to overburden the result, we restrict ourselves to 
a non-degenerate energy band, i.e. with one-dimensional 
eigenspaces. 

Theorem 2 In addition to Assumptions[^and[E assume 
that Ef is non-degenerate and that arbitrary derivatives 
of the corresponding family of eigenf unctions Lpi are glob- 
ally bounded. 
Up to terms of order e^ the second-order effective Hamil- 

(2) 

tonian H^f^ is given by 



{^' 



l^iffV^ 



{g^pfrpfr + vMri^ 



wher 

ffoff 



Pf 



-e^'ilj''U^RHiUiijndfi, 
+e^Wl,ig'"^Wi,^{cpf\3n'^n^^f), 



-isd, - e{ipi\iV)ipf) - e^TZ^JJiffl^n^n^id. 



'y^'Ji^f 



2 in" 



W 



(^f|n>f))iV^^f), 



Kff = Pf + e(9„l^)„=o(</'f|ri"v^f> 



-e' Kc 



VbH + Kmb + W2) 



Ui 



25'^ v^Vf^j + n"(ptW5'a?, - n'-ip{id^w)n=o 



and 

Vkmb 
W2 



2 I] 



Ti 



ab 



= <?^^"(v^Vf|(i-l^f)(vfl)v>f) 

= TZ\ p{d^Lpf\ln"n'^dsff}, 

= {dlf,W)n^o{flKn^cpf) 



7^ 



aj 



n- 



with W the Weingarten mapping, rj the mean curvature 
vector, TZ and TZ the Riemann tensors of C and A (see 
Appendix 1 for the definitions) . 

This effective Hamiltonian is derived in [4]. One might 
wonder whether the complicated form of the effective 
Hamiltonian renders the result useless for practical pur- 
poses. However, as explained in the introduction, the 
possibly much lower dimension of C compared to that of 
A outweighs the more complicated form of the Hamilto- 
nian. Moreover, the effective Hamiltonian is of a form 
that allows the use of semiclassical techniques for a fur- 
ther analysis. Finally, in practical applications typically 
only some of the terms appearing in the effective Hamil- 
tonian are relevant. As an example we discuss the case 



of a quantum wave guide in Section IV At this point we 
only add some general remarks concerning the numerous 
terms in H^f^ and their consequences. 

• The off-band coupling UlRHfUi can easily be 
checked to be gauge-invariant, i.e. not depending 
on the choice of Lpi but only on Pq. It occurs due 
to the replacement of Uq by U^ and thus is missed 
when one expands the naive adiabatic Hamiltonian 
i?^g — UoPoHePoUQ. Even if one, in addition, uses 
standard perturbation theory in the fibers, still the 
first term in Ui, which originates from [— e^Ah, Pq], 
would be missing. 

• Both T^BH , an analogue of the so-called Born-Huang 
potential, and Vamb, already found in P], are also 
easily checked to be gauge-invariant, which justifies 
to call them extra potentials. 

• The occurence of the geometric potential V^oom has 
been stressed in the literature, in particular as the 
origin of curvature-induced bound states in quan- 
tum wave guides (reviewed by Duclos and Exner 
in |26]), see Section [lV| below. 



E. Approximation of bound states up to order e 

(2) 
The unitary equivalence of H^ and -ff^g up to errors of 

order e^ allows us to deduce that the lower parts of their 

spectra coincide up to errors of order e^ when Ef is the 

ground state band. The following result, which is proved 

in [T, shows how to obtain quasimodes of iP^ from the 

(2) 
bound states of H^g and vice versa. 
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Theorem 3 Let Ef be a non-degenerate constraint en- 

(2) 

ergy band and let U^,H^g be the operators associated 
with Ef in the preceding subsection. 

a) Let E G M.. Then there is a C < oo such that for 
any family (Eg) withlunswp^^Q E,; < E and all e small 
enough the following implications hold: 






liHe- 



Es)u:^Pe\ 



<Ce'\\U:i:, 
II <Ce^W\\ 



b) Let E{{q) = inf (T(iJf (g)) for some (and thus for all) 
q E C and define Ei{q) := inf (cr(iJf (g)) \ i?f(g)). Let 
{ip^) be a family with 



limsup(?/''^|iff-0^) < infill. 



(30) 



Then there is c > such that WU^tp'^W > cH?/;'!! for all e 
small enough. 

We recall that for any self-adjoint operator H the bound 
\\{H - X)'ip\\ < 6\\iIj\\ for a e M implies that H has spec- 
trum in the interval [A — (5, A + 5] . So a) i) entails that H'^ 
has an eigenvalue in an interval of length 2Ce^ around E^, 
if one knows a priori that the spectrum of H"^ is dis- 
crete below the energy E. The statement b) ensures that 
a) ii) really yields a quasimode for normal energies below 
infqgc^'i: i-G- that 






Ee)U,^//\\ < ^e'llC/eV' 



If the ambient manifold A is flat, then (|30|) follows from 



limsup(^^|i/e^''> < inf £;i 



sup {Wn=0 

(q-.n) 



W) =: E, 



(31) 

Therefore Theorem [3] in particular, implies that at least 
for flat A there is a one-to-one correspondence between 
the spectra of H^ and H^g below £'*. In the example 



of Section I B depicted in Figure [2] this implies that all 
eigenvalues of H^ in the interval [eo,ei) and the corre- 
sponding eigenfunctions are determined by the effective 
Hamiltonian of the ground state band Eq modulo terms 
of order e^. 

(2) 

The bound states of H^f^ can be approximated by the 
standard WKB construction. In the simplest case one 
obtains: 

Corollary 1 Assume that A is flat and that Ef is a non- 
degenerate constraint energy band with inf _Ef < E^ and 
Ef{q) = inf cr(i/f (g)j for all q E C. Let there be qo E C 
such that Ef{qQ) < Ef{q) for all q y^ qo and (Vf £'f)(qo) 
is positive definite. 

Denote by Eg{A) the i-th eigenvalue of a semi-bounded 
operator A, counted from the bottom of the spectrum. 
Then for any £ G N 

Ee{H') = Efiqo) + sEeiHw) + ^(e'), 

where Hhq := — Ajjd -I- |(V?, . g .Ef){qo)x'^x^ is a har- 
monic oscillator on M.'^ . 



IV. QUANTUM WAVE GUIDES 

In this section we look at the special case of a curve C 
in A — M.^ equipped with the euclidean metric. Such 
curves may model quantum wave guides which have 
been discussed theoretically for long times (see e.g. the 
review pS]) but are nowadays also investigated experi- 
mentally (see e.g. the review [3?!). 

In the first subsection we provide the expression for our 
effective Hamiltonian when applied to wave guides and 
make some general remarks about trapping and splitting 
of wave packets. In the second subsection we explain how 
to produce topological phases in closed wave guides. The 
effects on the spectrum of such wave guides are discussed 
in the last subsection. 



A. Trapping and splitting in quantum wave guides 

We first look at infinite quantum wave guides. So let the 
curve C be given as a smooth injective c : K — >■ K'^ that 
is parametrized by arc length (|c| = 1). The mean curva- 
ture vector of c is ?7 = c and its (exterior) curvature is \r]\. 
Denoting by • the usual scalar product in K"^ we define 
y{n) := n ■ r]/\r]\ where ry 7^ and y{n) := elsewhere. 
By the Frenet formulas the Weingarten mapping satisfies 
y^iv) — hP (see e.g. |28]) and W = on the orthogonal 
complement of rj (which is NqC if rj(q) = 0). 
A normalized section of the tangent bundle TC is given 
by T := c. We extend this to an orthonormal frame of 
TC X NC, where NC is the normal bundle, in the following 
way: We fix g G C, choose an arbitrary orthonormal basis 
of NqC, and take vi, 1^2 to be the parallel transport of this 
basis with respect to the normal connection V"*" (defined 
in the first appendix) along the whole curve. This yields 
an orthonormal frame of NC. Together with t we obtain 
an orthonormal frame of TC x NC, which is sometimes 
called the Tang frame. We denote the coordinates with 
respect to r, i^i, and 1^2 by x, ni, and 722 respectively. 
In these coordinates it holds V*^ — dx (as can be seen 



from the coordinate formula (12 1 and the definition of 
the connection coefficients w in Appendix 1). 

Now let Ef and (</?/) j be as in Assumptionl2l We start by 
spelling out the formula for H^^ from Theorem 1 Since 
C is one-dimensional and contractible, the families of (pf 
can be chosen such that p|g = —iedx globally. Then the 
first-order effective Hamiltonian is 






-iedx{l+e\rM(Pi\y^-f))iedx + Ef 
-£(a„T4^)„=o(vfl""¥'/) (32) 



with {(f)\ip) := L2 </** "0 dnidn2. 

For highly oscillating states '0, i.e. with (?/'|— e^9^^0) ^ 1, 
the only term of order 1 besides —e^d^x is Ef. So if Ef 
is constant, in particular, if the wave guide has constant 
cross section, the dynamics is free at leading order and, 
even more, the potential terms are of order e^. So they 
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only become relevant for times of order £~^. However, a 
semiclassical wave packet ip covers distances of order e~^ 
on this time scale. Hence, for such tp noteworthy trapping 
occurs only for very long wave guides! 

If we consider a straight wave guide, i.e. f? = 0, the for- 
mula we end up with is the expected adiabatic approxi- 
mation: 



-'-'qwgl'7=0 



-292 



El + e{dc.W)n=o{'pi\n'- 



fi) 



We note that, although 77 = 0, the x-dependence of the 
constraining potential still allows us to model interesting 
situations. For example a beam splitter may be realized 
by fading a single-well into a double-well potential (see 

e.g. m)- 



B. Topological phases in quantum wave circuits 

Up to now we have considered an infinite wave guide, 
which, of course, has the topology of a line. The only 
possible non-trivial topology for a curve C is that of a 
circle. We refer to a wave guide modeled over such a C as 
a quantum wave circuit. In order to keep formulas sim- 
ple and transparent, we look at a so-called round circle, 
that is with constant 77. Then the Tang frame from the 
preceding subsection is still globally smooth. However, 
because of the non-trivial topology our choices of the 
families Lp'^ made above are only possible locally but in 
general not globally. Therefore we rewrite ( 32 1 without 
those choices. For the sake of brevity, we assume that W , 



the non-constraining part of the potential, is identically 
zero in the following. 



qwc 



pls{l + e\r^\{ipi\v^i))p,s + Ei (33) 



withpeff = — ie^a^+e ('/'f |i9a;<^/). Although the curvature 
of the connection ipeff always vanishes, it may lead to a 
topological phase, which we will discuss next. 

Here and in the following subsection we again restrict 
ourselves to the case of a non-degenerate energy band Ef. 
We note that even for degenerate energy bands only 
abelian phases will occur because the fundamental group 
of the circle is generated by only one element. Let x be a 
27r-periodic coordinate on the circle. The eigenfunction 
(/?f (x) associated to E^ can be chosen real- valued for each 
fixed X because Hi is real. This associates a real line bun- 
dle to Ef. From the topological point of view, there are 
exactly two real line bundles over the circle: the trivial 
one and the non-trivializable Mobius band. In the former 
case the global section ipi can be chosen real everywhere. 
This implies 

(vPfl^.V'f) - \{{Vf\d.^f) + (5.¥>f|^f)) = ^^^i^^ ^ 0, 

which results in ipcs — sdx- Thus there will be no topo- 
logical phase in this case. We will now provide an exam- 
ple for the realization of the Mobius band by a suitable 




FIG. 3; Level set of a potential like Vc describing a 
"Mobius wave circuit" . 



constraining potential and show that, indeed, a topolog- 
ical phase occurs! 

Let Vc G C^{M.'^) have two orthogonal axes of reflection 
symmetry, i.e. in suitable coordinates 



Vc{-ni,fi2) = Vcini,h2) = Vcini,~n2). 



(34) 



Then the real ground state $0 of — Aa2 + Vc with en- 
ergy Eq is symmetric with respect to both reflections, 



$0(^1,^2) == ^Q{-ni,n2) = $0(^1 



'"-2) 



while the first excited state $1, also taken real- valued, 
with energy Ei is typically only symmetric with respect 
to one reflection and anti-symmetric with respect to the 
other one, e.g. 

$1(^1,^2) = -*i(-ni,n2) = $1(^1,-^2). (35) 

This is true in particular for a harmonic oscillator with 
different frequencies. As the potential constraining to the 
circle we let Vc perform half a twist along the circle, i.e. 

{Vc^/^{x)){nun2) := I4(cos(x/2)ni ~ sin(x/2)n2, 

sin(x/2)r7,i -|- cos(x/2)n2). 



=/2 



We note that due to (34) this defines a smooth Vc 
Then 

((^j(x))(ni,n2) :— <f>j( cos(x/2)ni — sin(x/2)n2, 

sm{x/2)ni + cos(x/2)n2) 

is an eigenfunction of H{{x) := — Av -I- Vc{x) with eigen- 
value Ej for every x and j e {0, 1}. However, while 1^0 
is a smooth section of the corresponding cigenspace bun- 
dle, (fi is not. For by (35) it holds (pi(x) = — (^i(a; -I- 2tt) 
(see Figure [4]) . Still the complex eigcnspace bundle ad- 
mits a smooth non-vanishing sect ion. A possible choice 
is iy5i(x) := e'^/2|^i(x). Using (35) we obtain that for the 



first excited band the effective Hamiltonian ( 33 1 reduces 
to 

i?qwc,i = (-ie5, + e/2)2 + E^, 

while for the ground state band it is 



H, 



qwc,0 



2^2 

XX 



-e^d. 



Eo- 
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This shows that depending on the symmetry of the nor- 
mal eigenfunction the twist by an angle of tt has different 
effects on the effective momentum operator in the effec- 
tive Hamiltonian. With respect to the connection ap- 
pearing in iJq^j. I the holonomy of a closed loop 7 wind- 
ing around the circle once is h{-/) — e'Jo 1/2^3; _ _y 
Hence, the 1/2 cannot be gauged away. Furthermore, a 
wave packet which travels around the circuit once accu- 
mulates a topological phase equal to tt. 




x = 



FIG. 4: First excited transverse eigenfunction in a "Mobius 
wave circuit" : ni and n^ are the normal elements of a global 
Tang frame and with respect to this frame the confining po- 
tential Vc twists by an angle of n when going around the cir- 
cuit once. Due to its symmetry Vc is still globally smooth, 
see Figure [3] However, since the first excited state "I>i has 
not the full symmetry of the potential, it changes sign after a 
twist by TV. 



C. Effects of twisting and bending on the spectrum 

At second order in e the effect of the topological phase 
in quantum wave circuits can also be seen in the level 
spacing of ffqwcj and thus, with Theorem [3] also in the 
spectrum of H^. So we add the corrections of second 
order from Theorem [2] to (33 1. Of course, all terms con- 
taining the inner curvature of C and ^ = M'^ vanish due 
to the flatness of C and A with the euclidean metric. 



qwc 



P*cB{^ + £\v\{Vf\yff) +3£^|77p(v3t|yVf»Pcff 
+ ARe \r]\ edx {d^(pi\RHtyff) £^d. 



.2 02 



,2g2 
-^xx 



(36) 



with poff given by 

-isdx + £{^pi\\dx^i) -t- £^\'q\{tfi\2 (y - {ipf\yipf))dx<fii). 

If Vc does not change its shape but only twists, Ei is 
constant and thus may be removed by redefining zero 
energy. Furthermore, since the remaining potential terms 
are of order e^, the kinetic energy operator — e^c*^^ will 
also be of order e^ at the bottom of the spectrum. So 



i?q^^ may be devided by e^. Keeping only the leading 
order terms we arrive at 



HqwcJ —PctfPcff^ 



i'?r 



+ {dx(Pj\dx^j) - 1 {ipj \dxfj)\'^ (37) 



with pcs :~ — i9j; + (^Lpjlidx'fij)- A simple calculation 
yields 

{dxfjldx'Pj) - |(v5j|5:r(^j)P 

We note that the integral is the expectation value of the 
squared angular momentum of $f and thus vanishes for 
a rotation- invariant $f. So (37) shows that bending is 



attractive, while twisting is repulsive. 

Since \r]\ is constant, for £ e No the £-th eigenvalue of 



SK^qwca)=^i+£'[(^+5)' + 
while for i/qwc,o we find 

Ei{Hq^c,o) = Eq + 6 £ + 



LH-S'o)-\n\- 



+ Oie'), 



+ 0(e3). 



We note that, although a constraining potential that 
twists along a circle was investigated by Maraner in detail 
in [2] and by Mitchell in [3j, the effect discussed above 
was not found in both treatments. The reason for this 
is that they allowed only for whole rotations and not for 
half ones to avoid the non-smoothness of ipi . 

There is a wide literature on the spectrum of a quan- 
tum wave guide which is arbitrarily bent and twisted 
(see the review J3_0J by Krejcifik). In general, the twist- 
ing assumption means that there is e Cj^(M) and 
Vc e Cj^(M^) such that the constraining potential has 
the form: 

{Vc{x)){ni,n2) := V'c('^i cos6'(a;) - ^2 sin6'(x), 

ni sin d{x) + 712 cos 9{x)) . 

Then the family of eigenfunctions tpi may be chosen as 

((^f(a;)) (711,712) := ^f(niC0s6{x) — n2Sm9{x), 

rii sin 9{x) + 712 cos 0{x)) 

for an eigenfunction $f of — A^a + Vc{x) with eigen- 
value E[. It is easy to generalize the discussion above 
to a wave circuit whose curvature and potential twist are 
non-constant. Then the i?-th eigenvalue of H^ is given by 

EeiH,) = Et + e^i^K^twist) + 0{e'), 

where Eg{H^^-^^^) is the l-i\v eigenvalue of the following 
operator: 



Hi 



p:gPeff-hiV4+i'($f)^V4. 
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with poff := -iSx + {(pi\\{dx + {va ■ 9a;i^/3)""9^)</'f)- This 
generahzes results by Bouchitte, Mascarenhas and Tra- 
bucho [3T] and by Borisov and Cardone [32] for wave 
guides to wave circuits. We note that for the Tang frame 
Va'dxVp = but for an arbitrarily curved wave circuit the 
Tang frame is, in general, not globally smooth anymore. 
Anyway the normal bundle is still trivializable because it 
inherits the orientation of M^ and every orientable vector 
bundle over a curve is trivializable. 



V. CONCLUSIONS 

While all earlier results on constrained quantum systems 
had to focus either on a certain energy regime or on spe- 
cial geometries, we have presented here results, both on 
the dynamics and on the spectrum, that cover all relevant 
energy regimes in general geometries (recall Figure [2]) . 
We point out that our results on dynamics (Theorem [I] 
and Theorem [2]) are true for all bound state and scatter- 
ing energies, as long as oscillations faster than e~^ are 
excluded. The same is true for the quasimodes of the full 
Hamiltonian H'^ constructed from those of the effective 
Hamiltonian (Theorem Isl) . 

Furthermore, we have applied our results to quantum 
wave guides and obtained for the first time the complete 
second order effective Hamiltonian ([36]). In contrast to 
earlier theoretical results it applies also to wave circuits, 
i.e. closed wave guides. Here the effect of an abelian 
topological phase is observable both in the spectrum and 
in the dynamics. We believe that as a next step it would 
be interesting to apply our results to simple examples 
from molecular dynamics, like those that were treated for 
small kinetic energies by Maraner in [2T]. Here also the 
curvature of the effective Berry connection, calculated in 
Proposition [T] should play a role. Note that it did not 
show up for quantum wave guides because of the one- 
dimensional constraint manifold. 



Appendix 

Al Geometry of submanifolds 

We recall here some standard concepts from Riemannian 
geometry. For further information see e.g. [33]. As be- 
fore we use the abstract index formalism including the 
convention that one sums over repeated indices. More- 
over, we will consistently use latin indices a,b, .. running 
from 1 to d -|- fc for coordinates on a general manifold, 
latin indices i, j, .. running from 1 to d for coordinates on 
a submanifold, and greek indices a, /3, . . . running from 
d -|- 1 to d -|- fc for coordinates in the normal spaces of a 
submanifold. 

First we give the definition of the Riemann tensor we use 
(the order of the indices varies in the literature!). 



Definition 1 Let {A4, g) be a Riemannian manifold with 
Levi-Civita connection V. Let (tq)^ be a set of local co- 
ordinate vector fields, 
i) The Christoffcl symbols F 0/ V are defined by 

ii) The Riemann tensor TZ is given by 

T^lcd := (VcVdTb-VdVcTb)". 

As usual by raising and lowering indices we mean to shift 
covariant to contravariant coordinates and vice versa. 

Now we turn to the basic objects related to the exterior 
curvature of a submanifold of arbitrary codimension. 

Definition 2 Let (A, g) be a Riemannian manifold with 
Levi-Civita connection V. Let C <Z A be a submanifold 
equipped with the induced metric g = g\c. Denote by NC 
the normal bundle of C . Let {Ti)i be a set of local coor- 
dinate vector fields of C and {va)a o, local orthonormal 
frame of NC . 
i) The Weingartcn mapping W is given by 

ii) The second fundamental form II is defined by 
Hi) The mean curvature normal r/ is defined by 



Va 



Wi 



The relations and symmetry properties of W and II for 
hypersurfaces also hold when the codimension is greater 
than one: 



IP . = WV 



W: 



Hi 



(38) 



Finally, we provide the definitions of the objects that 
characterize the geometry of the normal bundle. 

Definition 3 i) We define the normal connection V^ to 
be the bundle connection on the normal bundle given via 

for a = d+1, . . . , d-{-k and {Wj-Ua)"" = for a — 1, . . . ,d. 
ii) The connection coefhcients lu of V^ are defined by 

Hi) The normal curvature tensor R^ is defined by 

Due to the anti-symmetries of any curvature tensor the 
normal curvature tensor R^ is identically zero, when the 
dimension or the codimension of C is smaller than 2. 

When we set W„" = for a = d+ l,...,d + k, the 
Weingarten equation 

[vfv^r = {v,u^r+w^^^ (39) 

is a direct consequence of the definitions. 
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A2 Expansion of the metric tensor 

In order to expand the Hamiltonian H^ in powers of e 
it is crucial to expand the metric g on the normal bun- 
dle NC around C because the Laplace-Beltrami operator 
depends on it. Such expansions were carried out in al- 
most any work on constrained quantum systems, how- 
ever, in various generalities and up to varying orders. 
Here we provide a simple derivation for an arbitrary sub- 
manifold of a curved ambient manifold but only up to 
first order. However, this is enough in order to obtain 
Theorem [TJ 

Fix (J G C Let {Ti)i be a set of local coordinate vector 
fields of C and {va)a a local orthonormal frame of NC. 
Furthermore, let ^(q) := {(f)^{q), . . . , (j)'^^'^{q)) be the lo- 
cal expression for the isometric embedding of C into NC 
and $((7, N) be the exponential map in each fiber. Then, 
by definition of the exponential map, $((7, N) = 7(1) 
where 7(5) is the geodesic on NC starting at <j){q) with 
7(0) = N'^i^^. 

Let (y°)a=i,...,d+fe be Riemannian normal coordinates on 
{NC,g) around (f){q) and FJJ^ the associated Christof- 
fel symbols of the Levi-Civita connection. Due to the 
geodesic equation 7°(s) = —F^^ (7(3)) 7'' (5)7^(5) a Tay- 
lor expansion around s = yields 

7"(s) = riq) + sN"i^^iq) 

-^N'^N0nMq)y^{q)M^^iq) + 0{s^) 

Evaluating at s = 1 we obtain that 

-i7V"iV^F^,(0(g))^^(9)z.^(g) + 0{\N\'). 

Therefore 

d,<P''{q,N) = d,r{q) + N"d,i^l{q) + 0{\Nf) 
d^^'^iq,N) =. «-i7V^F?,(0(.))^^^)((?) + O(|iVn 

= <(g) + 0(|iVp), 

where we used that F^^((/)((7)j = in Riemannian 
normal coordinates. The latter also implies that 
C^ i'^a)"' {4>iq)) = diiy^{q). Then the Weingarten equa- 
tion ( 39 ) yields that 



d^Mq) = Loliq^^iq) - WUq)dl^iq)- 

Using that 9,0 and Wj^idicj) are tangent vectors and thus 
orthogonal to vp for any i and /3 we obtain that 

9^j{q,N) = {g,,d,<^''d,<i>'){q,N) 

~9abN"W'a^^lr^,^')iq)+Oi\N\') 
g,p{q,N) = {g,,d,<^''dp<^''){q,N) 

= {Va.N'^^l.u^A) (9) + 0{\N\^) = g^,(g, N) 
g„p{q,N) = {g,,d^<^'^dp<^''){q,N) 

= {9a,<4)^q) + 0{\N\^), 



Since is an isometric embedding, it holds 
g^i,di(j)'^dj4>^ = gij. The orthonormality of the nor- 
mal frame yields Vat'^a'^B ^ ^ap- Thus 

g,^iq,N) ® g,,.(g)-2iV"H„,,.(g)-HO(|A^n, 
g,piq,n) = 7V"c.,„^(g) + 0(|7V|2) = gp,iq,N), 
g^0{q,N) = 6^p + Oi\N\^). 

Inverting this matrix we end up with this proposition: 

Proposition 3 The inverse metric tensor g has the 
following form for all q G C: 



g{q,en) = 



1 
C^ 1 



A 
B 



1 ) (9>e")> 



where for i, j,l,m ~ 1, ...,d and a,/3,^,6 — d+l,..^d + k 
A^\q,en) = g^^iq) +s2n'^lV^{q) + 0{eM^), 

CUq,en) = ^sn'^ uUq) + 0(6^')- 

Here II is the second fundamental form and lo are the 
coefficients of the connection on the normal bundle (see 
Appendix 1 for the definitions). 

We note that the error is of order e^ jrip not only for small 
\n\ but globally, when the metric g is chosen properly 
outside a tubular neighborhood of C (see Section |II A 
and |4]). Then the use of the expansion is justified by 
the fast decay of functions from the subspaces Pq and P^ 
in the fibers. 



A3 Transformation of measures 

Let (Ti be the density of the measure djl on NC and (T2 
be the density of the product of the measure d/x on C 
and the Lebesgue measure dN on the fibers NqC = M'^. 
Define p := — and 

r cr2 

Mp : L^{NC,d-p) ^ L^{NC,dfi(g) dN), ^ ^ p-^. 
Mp is an isomctry because for all *, # e L^{NC, dfKi^dN) 



Mp^Mp^dJI^ / -^-^ p^'^d-p= I'l'd^^dA^. 

NC JnC JnC 



Therefore it is clear that 
One immediately concludes 



MpM*p = 1 = m;Mp 
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and thus Mp is unitary. Now we note that [db^p 2] = 
— \ /9~2 9f,lnp. So for AdcTi '■— -'(^i'^da(Tig°'^db we have 

M;(-AdajA/p 

On the one hand, 

p^-a:['^da(Tip-^ g^^db 

= pa^'daaip-^g'^'db + y''\da\np)db 
= <J2^daa29'"'db + y''''{db\np)da 

and on the other hand, 

p^a^^daaip^^g°-''dblnp 
= ~y'^'{dalnp)idb\np) 

ma^'daaig^^'dblnp) + y''''{db\np)da. 
Together we obtain 
M;(-Arf,jMp 
= -<y2'da<J2g''''db 

~ y''\da In p) (db In p) + i {a^'daCJig''''db In p) 
= -Ada.i^ - \g''''{da\np){db\np) + U^dajnp). 



Because of Aj^^--^ — A^ we have shown that 



M;{-Ag)Mp 



-Ada, + Vp 



(40) 



with Vp := -lg''''idalnp){dblnp) + i(Ad,,hip). This 
formula was estabhshed many times before and we have 
provided its derivation for the sake of completeness, as it 
is the origin of the geometric potential. 



A4 Expansion of the Hamiltonian 



In order to deduce the formula for the effective Hamilto- 
nian we need that H^ = — e^(Ag)'^ + Vdq, n) + W{q, en) 
can be expanded with respect to the normal directions 
when operating on functions that decay fast enough. For 
this purpose we split up the integral over NC into an in- 
tegral over the fibers NqC, isomorphic to M*^, followed by 
an integration over C. 
The following expansion is also the justification for the 



splitting of H^ in ( 16 1 



Proposition 4 // an operator A satisfies 

\\A{nY\\ciu) < Ci, \\{nyA\\c(v(H^)) < Q 

for every I G N, then the operators H^A and AH;, can be 
expanded in powers of e on C(T>{H^),H): 

H,A ^ {Ho + eHi)A + 0{e^), 
AH, = A(Ho + eHi) + 0{e^), 



where Hq and Hi are the operators associated with 

Jc Jr'' 
(*|ffi*) = f f 2n"iri(eVh*)eV^'* (41) 

Jc JR'' 

-fn"(a„l¥)„=o|*Pdrid/i, 

where \/^ is the horizontal connection (see Section 
for the definition). 



HB 



To derive this let P with \\{n)'P\\c(v(H,)) < Ci for I e Nq 



be given. The similar case of a P with \\P{ny\\ciH) ^ ^i 



for all I e No will be omitted. 

We set *p := P*. By definition of H, it holds 

(*|i/e*p) = (*|-e'(Ag)^*p) 

+ (*| (K((7, n) + W{q, en))^p). (42) 



The formula ( 40 ) implies that 



(* I - e^Ag^p) 



e^ g''''da-9db'^ P + e^Vp^-^ P dN dp 



e' g'^'da'^db-^p dN dp + ©(e^). 



We emphasize once again that the remaining term may 
be of order 1 for a ^ with energy of order 1. To calcu- 
late — e^(Ag)'^ we have to replace N by en in the above 
formula. Then we may exploit \\{n)'^ P\\c{v{h^) ^ C2 to 
insert the expansion for g from Proposition p] into ( 43 1 . 
Noting that the rescaling n — N/e does not change9i 
and replaces da by £"^9q we obtain that 



{^\-e^A^r^p) 



c Jr'' 



{ed, + Cf{q,en)da)'i>A'^{q,en) 



x{ed,+C^^iq,en)dppp 
+ ~{da^B"^{q,en)dp^pdndp + 0(6^) 

(£(9,-n7wf^(q)5„)*) {g''+e2nnV^ 

xe(a,-nVf^((7)a^)*P 
+ Jda^6°''^dp^pdndp + 0{e^) 



(eV^'*) {g'^ +e2n'^lVi)ey)-^p 

+ *"(-A„)*pdrid^ + 0(£2),(43) 



where we used (12 1 and A„ — S'^^dadn. Due to 



||(n) P|| < C2 a Taylor expansion of W{q,en) in the 
fiber yields 

W{q, en)P = {W{q, 0) -f- en"{daW){q, Q))P + 0{e'^). 



Plugging this and ( |43| into (W2| we obtain the claim when 
we recall the definition Hf :— — A„ -I- Vc{q, n) -{- W{q, 0). 
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A5 Derivation of the effective Hamiltonian 



with 



Here we derive the formula for H^f^ ~ UqH^Uq stated in 
Theorem [T] Plugging in the expansion of i/^ from the 
preceding appendix we have that 



„JK 



H. 



(1) _ 






In the following we write {■ \ -Jnc for the scalar product 
on L'^{NC,dNdfi). By Definition of Uq in Theorem [l] we 
have 

{r\UoAU*ij')c = {y^iri\A^Uj)NC- (44) 

for any oper ator A. In view of the definition of V in 

Section [lIB 

connections: 



V*^ satisfies the usual product formula for 



eVl'^^^f = ^'fEd.ri+eri^'^^l (45) 

We note that eV^t^f is really of order e, while edi^'j is, 
in general, of order 1 due to the possibly fast oscillations 
of '0'^. Furthermore, the exponential decay of the (pf, 
which implies the exponential decay of their derivatives 
(see [?j), guarantees that, in the following, all the fiber 
integrals are bounded in spite of the terms growing poly- 
nomially in n. The product formula (45) implies that 



{ipjipj\Hoipfip';j)NC 









u: 



g^'SuPi^.^KPi^j^l + EiS^'^^j^jdndii + 0{e 



p;^-^ = -ie<5^^a,-£(^/|iV>f). 



Furthermore, 

{ifi^pj I Hi ipWj)NC 

JIT 



/ / 2n"II^J(eVh(^fVf)eV^</5/V'jdnd/z 
Jc Jv' 



n°'{daW)„=o(pj'iljliPf^pjdnd^i 



2nnV^ {(fijed.i^^j) t/s/eSj V'j dn dfi 



+ / {daW)n=o{v'fK^'l) ri^'j dfi + 0{e) 



21V^ {vi\n"^i)pi^rKPiii^ 



li'FL 



+{dc.W)n=o{^'i\n"ipf)i;frj dfi + 0{e), 
where we used that p;;!.^, = —ieS'^^dj + 0{e). So we, 
indeed, obtain that 



xU 






with 



g%" = g^^S" + elV^{^Un'^^^). 
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